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Motivated by the recent work of one of us [lj], we generalize this work to the case where 
the pressureless dark matter and the holographic dark energy do not conserve separately 
but interact with each other. We investigate the cosmological applications of interacting 
holographic dark energy in Brans-Dicke theory with chameleon scalar field which is non- 
minimally coupled to the matter field. We find out that in this model the phantom crossing 
can be constructed if the model parameters are chosen suitably. We also perform the study for 
the new agcgraphic dark energy model and calculate some relevant cosmological parameters 
and their evolution. 
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I. INTRODUCTION 

Among various scenarios to explain the acceleration of the universe expansion, the holographic 
dark energy (HDE) and agegraphic dark energy (ADE) models have got a lot of enthusiasm recently. 
These models are originated from some considerations of the features of the quantum theory of 
gravity. That is to say, the HDE and ADE models possess some significant features of quantum 
gravity. Although a complete theory of quantum gravity has not established yet today, we still 
can make some attempts to investigate the nature of dark energy according to some principles of 
quantum gravity. The former is motivated from the holographic principle [3, [3] • It was shown in 
4{| that in quantum field theory, the UV cutoff A should be related to the IR cutoff L due to limit 
set by forming a black hole. If pr> = A 4 is the vacuum energy density caused by UV cutoff, the 
total energy of size L should not exceed the mass of the system-size black hole: 

E D < E BH -> L 3 p D < m 2 p L. (1) 

If the largest cutoff L is taken for saturating this inequality, we get the energy density of HDE as 

3c 2 m2 3 C 2 

PD ~ ~ir ~ ugi?- (2) 

The HDE is thoroughly investigated in the literature in various ways (see e.g [jj and references 
therein ). The later (ADE) model assumes that the observed dark energy comes from the space- 
time and matter field fluctuations in the universe. Following the line of quantum fluctuations of 
spacetime, Karolyhazy et al. [{J] discussed that the distance t in Minkowski spacetime cannot be 
known to a better accuracy than 8t = where j3 is a dimensionless constant of order unity. 

Based on Karolyhazy relation, Sasakura [7|] discussed that the energy density of metric fluctuations 



of the Minkowski spacetime is given by (see also 



PD ~ ^ ~ -p, (3) 

where t p is the reduced Planck time and t is a proper time scale. On these basis, Cai [s| proposed 
the energy density of the original ADE in the form 

PD = —^2*, (4) 

where T is the age of the universe. Since the original ADE model suffers from the difficulty to 
describe the matter-dominated epoch, the new ADE (NADE) model was proposed by Wei and Cai 
[lo| . while the time scale was chosen to be the conformal time instead of the age of the universe. 
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The APE models have arisen a lot of enthusiasm recently and have examined and studied in ample 
detail 
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It is also of great interest to analyze these models in the framework of Brans-Dicke (BD) gravity. 
In recent years the BD theory of gravity got a new impetus as it arises naturally as the low energy 
limit of many theories of quantum gravity such as superstring theory or Kaluza-Klein theory. The 
motivation for studying these models in the BD theory comes from the fact that both HDE and 
ADE models belong to a dynamical cosmological constant, therefore we need a dynamical frame to 
accommodate them instead of Einstein gravity. The investigation on the HDE and ADE models in 



15H18I] . In the present work, we consider 



the framework of BD cosmology, have been carried out in 
a BD theory in which there is a non-minimal coupling between the scalar field and the matter field. 
Thus the action and the field equations are modified due to the coupling of the scalar field with 
the matter. This kind of scalar field usually called "chameleon" field in the literature [3]. This 
is due to the fact that the physical properties of the field, such as its mass, depend sensitively on 
the environment. Moreover, in regions of high density, the chameleon blends with its environment 



and becomes essentially invisi 



)le to searches for Equivalence Principle violation and fifth force [jjj] . 



Further more, it was shown Ijjj [^(J that all existing constraints from planetary orbits, such as 
those from lunar laser ranging, are easily satisfied in the presence of chameleon field. The reason 
is that the chameleon-mediated force between two large objects, such as the Earth and the Sun, is 
much weaker than one would naively expect. In particular, it was shown [20( that the deviations 



from Newtonian gravity due to the chameleon field of the Earth are suppressed by nine orders of 
magnitude by the thin-shell effect. Other studies on the chameleon gravity have been carried out 



m [21 1 . Our work differs from that of Ref. [1711 in that we assume a non-minimal coupling between 

n 

the scalar field and the matter field. It also differs from that of Ref. [1[, in that we assume the 
pressureless dark matter and dark energy do not conserve separately but interact with each other, 
while the author of []J assumes that the dark components do not interact with each other. 

II. HDE IN BD THEORY WITH CHAMELEON SCALAR FIELD 



We begin with the BD chameleon theory in which the scalar field is coupled non-minimally to 
the matter field via the action 
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S = / tfxy[=g <j>R - -g^d^d^ - V{4>) + f{<f>)L 



(5) 
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where R is the Ricci scalar curvature, <j) is the BD scalar field with a potential V((f>). The 
chameleon field (f) is non-minimally coupled to gravity, u is the dimensionless BD parameter. The 
last term in the action indicates the interaction between the matter Lagrangian L m and some 
arbitrary function /(</>) of the BD scalar field. In the limiting case /(</>) = 1, we obtain the 
standard BD theory. 

The gravitational field equations derived from the action ([5]) with respect to the metric is 

where T^ v represents the stress-energy tensor for the fluid filling the spacetime which is represented 
by the perfect fluid 

Tfw = {p + p)u^u u + pg^ u , (7) 

where p and p are the energy density and pressure of the perfect fluid which we assume to be 
a mixture of matter and dark energy. Also is the four-vector velocity of the fluid satisfying 
u^Un = — 1. The Klein-Gordon equation (or the wave equation) for the scalar field is 

where T is the trace of (|7|). The homogeneous and isotropic Friedmann- Robertson- Walker (FRW) 
universe is described by the metric 



ds* = -dt 2 + a'(t) I ^frWI , (9) 

where a(t) is the scale factor, and k = — 1,0, +1 corresponds to open, flat, and closed universes, 
respectively. Variation of action (0) with respect to metric @ for a universe filled with dust and 
HDE yields the following field equations 



k u<j> 2 | H <f> _ f{4>) f _ , , v(<f>) 

a 2 6 (j) 2 4> 30 



H +^2 + H li = (pm + Pd)+-£t, (10) 



2^ + H 2 + * + ^ + 2H^ + ± = (11) 

a a 2 2(j) 2 <f> (j) (j) 2<j) ' 



where H = a/a is the Hubble parameter, pd, Pd and pu are, respectively, the dark energy density, 
dark energy pressure and energy density of dust (dark matter). Here, a dot indicates differentiation 
with respect to the cosmic time t. The dynamical equation for the scalar field is 

* + 3H * ~ W&(f ~ 5«*) + 5^5 (" - t> v *) = a (12 » 
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We assume the HDE in the chameleon BD theory has the following form 

The motivation idea for taking the energy density of HDE in BD theory in the form (|13p comes 
from the fact that in BD theory we have <p oc G _1 . Here the constant 3c 2 is introduced for later 
convenience and the radius L is defined as 

L = ar(t), (14) 

where the function r(t) can be obtained from the following relation 

■r(t) dr = r™<tt = Rh 
y/1 - kr 2 Jo a a ' 

It is important to note that in the non-flat universe the characteristic length which plays the role 

of the IR-cutoff is the radius L of the event horizon measured on the sphere of the horizon and 

not the radial size Rh of the horizon. Solving Eq. (|16p for the general case of the non-flat FRW 

universe, we get 

r(t) = -y=siny, (16) 

where y = y/kRh/a. Now we define the critical energy density, p CT , and the energy density of the 
curvature, pk, as 

p CT = 3(pH 2 , p k = (17) 
As usual, the fractional energy densities are defined as 



Q _ PM _ Pm ns x 
~ 3W (18) 



Per HH 

For latter convenience we rewrite Eq. (|20p in the form 



«* = — = 77^' ( 19 ) 

^ = - = 77^- (20) 



HL = ^=. (21) 



Taking derivative with respect to the cosmic time t from Eq. (j!4j) and using Eqs. (|16p and (|2ip 
we obtain 

L = HL + ar(t) = —7= — cos y. (22) 



Consider the FRW universe filled with dark energy and pressureless matter which evolves according 
to their conservation laws 



PD + 3Hp D (l + w D ) = 0, 
p M + 3Hp M = 0, 



(23) 
(24) 



where wd is the equation of state parameter of dark energy. At this point our system of equations 
is not closed and we still have freedom to choose one. We shall assume that BD field can be 
described as a power law of the scale factor, <fi oc a a . In principle there is no compelling reason for 



this choice. However, it has been shown that for small a it leads to consistent results 



15j |. A case 



of particular interest is that when a is small whereas uj is high so that the product aoo results of 



order unity 
bound on u) 



. This is interesting because local astronomical experiments set a very high lower 
23j; in particular, the Cassini experiment implies that ui > 10 4 



24. 



25J. Taking the 



derivative with respect to time of relation <fi oc a a we get 

4> = aHcf), (25) 
4> = a 2 H 2 (p + acpH. (26) 

Taking the derivative of Eq. (|13p with respect to time and using Eqs. (|22p and (|25[) we reach 

p D = H PD (a- 2 + 2^^ cos y] . (27) 



Substituting this equation in Eq. (J23J), we obtain the equation of state parameter 

w D 



1 2^ 

■-(a + 1) cosy. 

3 3c 



(28) 



It is important to note that in the limiting case a = (lj — > oo), the Brans-Dicke scalar field 
becomes trivial and Eq. (|28[) reduces to its respective expression in Einstein gravity 

1 2y/HE 



w D 



cos y. 



(29) 



3 3c 

We will see that the combination of the Brans-Dicke field and HDE brings rich physics. For a > 0, 
wd is bounded from below by 



l 

w D = --(a + 1) — . 

3 3c 



1 



(30) 

, the lower bound becomes 



Assuming Q,d = 0.73 for the present time and choosing c 
wd = — § — 0.9. Thus for a > 0.3 we have wd < — 1- This implies that the phantom crossing can 



be constructed in the BD framework. We can also obtain the deceleration parameter 



aH 2 



H 

w 



(31) 
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which combined with the Hubble parameter and the dimensionless density parameters form a set 
of useful parameters for the description of the astrophysical observations. Dividing Eq. (Ilip by 
H 2 , and using Eqs. (fT3l) . (l2Tj). p5j) and (l26j). we find 



1 



a + 2 



(a + l) 2 + a 



/aw 
V~2~ 



1 +Q k + 3{l D w D --Q\ 



(32) 



where the last term can be understood as a contribution of the potential energy in the total energy 
density i.e. 

V 



Substituting wd from Eq. ([28]) in (f32|) . we get 



Pc 



1 



a + 2 



(a + If + a - 1 ) + fc - (a + l)n D - cos y-£JV 



(33) 



(34) 



If we take CId = 0.73 and ~ 0.01 for the present time and choosing c = 1, aoj ~ 1, u; = 10 4 
and cosy ~ 1, we obtain g = —0.48 for the present value of the deceleration parameter which is in 
good agreement with recent observational results j^ . 

From equation (|12p . we can also estimate the mass of the chameleon field. This can be done by 
calculating the second derivative of the potential function with respect to scalar field 



1 



p-3p 



Following previous studies 
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281 ]. we choose 



V ^ = —^> ftt) = foe bo *. 



28J. We get 
(35) 

(36) 



Here M, /q and 6q are finite parameters whose values are model dependent. Making use of Eq. 
([36]) in ([35]h we obtain 

1 r M 4+n 6 o /oe 6 ° 



n- 



+ 



-(p-3p)(l 



(37) 



/)«+! ' 2 

Clearly when n — > (which corresponds to a constant potential), the mass of the scalar field will 
be dependent on the properties of f(4>). Moreover if only (f> = l/6o, the mass is determined by the 
scalar potential function alone. 



III. INTERACTING HDE IN BD THEORY WITH CHAMELEON SCALAR FIELD 

In this section we would like to construct a cosmological model based on the BD chameleon 
field theory of gravity and on the assumption that the dark energy and dark matter do not con- 
serve separately but interact with each other. Taking the interaction into account the continuity 
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equations becomes 



p D + 3Hp D (l + w D ) 
PM + 3Hp M = Q, 



■Q, 



(38) 
(39) 



where Q is an interaction term which can be an arbitrary function of cosmological parameters 
like the Hubble parameter and energy densities. The dynamics of interacting dark energy models 
with different interaction terms have been investigated in [291 ]. It should be noted that the ideal 



interaction term must be motivated from the theory of quantum gravity. In the absence of such 
a theory, we rely on pure dimensional basis for choosing an interaction Q. Hence following [^(J, 
we assume Q = Tpi> with T = 36 2 (1 + r)H where r = pm/pd is the ratio of energy densities and 
b 2 is a coupling constant. Note that T > shows that there is an energy transfer from the dark 
energy to dark matter. Combining Eqs. (|17p and (|25|) with the first Friedmann equation (jlOp . we 
can rewrite this equation as 



i + n k = f(4>)(n M + Vd) + n^+-n v , 



where 



Combining Eqs. ((2 
interacting HDE 



a 



/aui 



l 



(40) 



(41) 



1 



WD 



5 <« + i 



0J) and (|4ip with Eq. (|38[) we find the equation of state parameter of the 

b 2 



2JU 



D 



3c 



cos y 



i + n k + a ( 1 - — 



(42) 



In the absence of the BD field (a = 0, f{<j) 



1, V(4>) = 0), Eq. (|12|) restores its respective 



expression in non-flat standard cosmology 

1 2^/n D ~ 

w D 



■ cos y 



:i + ^fc). 



(43) 



3 3c "tip 

Next, we examine the deceleration parameter, q = —d/(aH 2 ). Substituting wd from Eq. (|42|) in 
Eq. ([32]) . one can easily show 



1 



a + 2 
3b 2 



/aw 



3/2 



(a + I) 2 + a - 1 ) + - (a + l)n D - cos y--O y 



i + n k + a [i- — 



(44) 



Comparing Eq. (144|) with (|34f) shows that in the presence of interaction the chameleon function 
f(4>) enters explicitly in q expression. This is in contrast to the usual BD theory where q of the 
interacting HDE model does not depend on the scalar field [l7]. 
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Finally we present the equation of motion of the dark energy. Taking the derivative of Eq. (|20p 
and using Eq. (j22H and relation = HQ' D , we find 

0^ = 20^-^-1 + ^008^, (45) 

where the dot is the derivative with respect to time and the prime denotes the derivative with 
respect to x = In a. Using relation q = — 1 — jp, we have 

Q' D = 2n D (q+— ^cosyV (46) 



where q is given by Eq. (|44p . This equation describes the evolution behavior of the interacting 
HDE in BD cosmology with chameleon field. 



IV. INTERACTING NADE WITH CHAMELEON SCALAR FIELD 



The above study can also be performed for the new agegraphic dark energy (NADE) model. In 
NADE, the infrared cut-off is the conformal time which is defined as 

pa da 



dt 

V= I — 

a 



(47) 



/o Ha*' 

In the framework of BD chameleon scalar field, we assume the following form for the energy density 
of the NADE 

3n 2 (b 



PD 



(48) 



where the numerical factor 3re 2 is introduced to parameterize some uncertainties, such as the 
species of quantum fields in the universe, the effect of curved space-time and so on. The respective 
fractional energy densities can be written as 

,2 



n D = p -R 



rr 



p cr H 2 7] 2 ' 

Differentiating Eq. (|48p and using Eqs. (|25p and (|49p we have 



Pd = Hp D [a 

na 



(49) 



(50) 



Substituting this relation in Eq. (I38p and using relations (I40p and (I4ip . we obtain the equation of 
state parameter of the interacting NADE 



w D 



1 2 

■1 a H 

3 3na 



i + n k + a 1 



~6~ 



1 



(51) 
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When a = 0, / = 1 and V = 0, the BD scalar field becomes trivial and Eq. (|5ip reduces to its 
respective expression in NADE in Einstein gravity [ijj]. From Eq. (|51|) . we see that even in the 
absence of interaction (b = 0), the the phantom crossing will take place in the the framework of 
BD theory provided the model parameters are chosen suitably. Indeed in this case (b = 0), wd can 
cross the phantom divide provided naa > 2\/TTd. If we take £Id = 0.73 and a = 1 for the present 
time, the phantom-like equation of state can be accounted if na > 1.7. For instance, for n = 4 and 
a = 0.5, we get wd = —1.02. When the interaction is taken into account the phantom crossing for 
W£> can be more easily achieved for than when resort to the Einstein field equations is made. 

In the context of BD chameleon scalar field the deceleration parameter of interacting NADE is 
obtained as 



1 

q 



a + 2 



( Q + 1)2 + a f™ _ A + ft fc _ ( Q + 3) o D + A^/ 2 _^ y 
V 2 / na 2 



i + n fc + a(i- — )--n v 



f{<l>) V V 6/2 

While the equation of motion for takes the form 



(52) 



n' D = 2n D ii + q-y—tL . (53) 

na 



V. CONCLUSIONS 



In this paper, we have considered interacting HDE model in the framework of BD cosmology 
where the HDE density pr> = 3c 2 /(87rGL 2 ) is replaced with = 3c 2 (j)/L 2 . With this replacement 
in BD theory, we found that the cosmic acceleration will be more easily achieved for than when 
the standard HDE is considered. Following the work of [l|, we assumed that the scalar field is 
non-minimally coupled with the matter field via an arbitrary coupling function f{4>)- In principle, 
the coupling between BD scalar field and matter field should be derived from a theory of quantum 
gravity. In the absence of such a theory, we have kept our analysis general regardless of the 



specification of f(cp). In the present paper, we have extended the work [l] by incorporating the 
interaction term in the HDE model. An interesting consequence of the present model is that it 
allows the phantom crossing of the equation of state of dark energy due to the presence of several 
free parameters. We have also performed the analysis for the NADE model and calculate some 
relevant cosmological parameters such as the equation of state, deceleration parameter and energy 
density parameter. 
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